Lowering of surface melting temperature in atomic clusters 
with a nearly closed shell structure 
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We investigate the interplay of particle number, N, and structural properties of selected clusters 
with 7V=12 up to N—562 by employing Gupta potentials parameterized for Aluminum and extensive 
Monte-Carlo simulations. Our analysis focuses on closed shell structures with extra atoms. The 
latter can put the cluster under a significant stress and we argue that typically such a strained 
system exhibits a reduced energy barrier for (surface) diffusion of cluster atoms. Consequently, 
also its surface melting temperature, Tg, is reduced, so that TJg separates from and actually falls 
well below the bulk value. The proposed mechanism may be responsible for the suppression of the 
surface melting temperature observed in a recent experiments. 

PACS numbers: 61.46.+W, 65.80.+n 
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I. INTRODUCTION 

The properties of small metal clusters have enjoyed 
considerable interest in recent years. Cluster properties 
can differ enormously from those of the bulk material 
due to the large surface-to-volume ratio and due to a 
remarkable sensitivity of the electronic structure to the 
cluster size and geometry.! These properties are of in- 
terest in technological applications, e. g. for catalysis. 2 
From a conceptual point of view clusters pose fundamen- 
tal questions in statistical mechanics of finite systems^ 

The melting process of small clusters is a complex phe- 
nomenon, the detailed rules of which are poorly under- 
stood. From early on it has been associated with iso- 
mer fluctuations 4 -^. More recent investigations of isolated 
Niia-xALc alloy clusters^ elucidate the relation between 
isomer fluctuations and the increase of entropy across the 
melting transition. A detailed overview about how the 
phase space grows with increasing particle number and 
the classification of isomers in terms of potential energy 
surfaces is given in Ref. 7j. A general overview of struc- 
tural properties of nanoclusters is found in Ref. [8|. 

By contrast, the understanding of the melting of very 
large clusters and bulk materials is considerably better 
developed. It is believed to be a strongly inhomogeneous 
process. In large, free metal clusters melting starts in 
the outermost atom shell, at the interface to the vac- 
uum, because there the thermal fluctuations of the par- 
ticle density have the lowest energy cost^iIS By feeding 
more heat into the cluster melting peals off shells from 
the solid cluster core, layer by layer. Following this logic, 
one concludes that melting in principle is a continuous 
transition - to the extent that each one of the individual 
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layers carries its own melting temperature. For atomic 
clusters this is the rule rather than the exception, because 
typically each layer has its own atomistic structure. This 
behavior is reflected in the temperature dependence of 
the specific heat, C(T), which is not necessarily a very 
sharply peaked function of temperature (as is the case for 
bulk samples), but rather may exhibit a strong inhomo- 
geneous broadening which reflects intcr-shcll averaging. 
A historical overview on "continuous melting" is given in 
Ref. HO. 

Depending on the crystal orientation^ 2 -, the surface 
may melt already at temperatures well below the bulk 
melting points. The thickness of the molten layer is 
strongly temperature dependent. It increases continu- 
ously with increasing temperature and eventually it di- 
verges - by definition - at the bulk melting temperature. 
(It is assumed here, that the thermodynamic state of 
the cluster interior is unique, i. e. it becomes indepen- 
dent of geometry in the thermodynamic limit.) These 
observations can be treated theoretically in model calcu- 
lations using different effective potentials 14 ' 15 as well as 
phcnomenologically-i£. 

Coming back to small systems, this dependency of 
melting on surface crystallography suggests that surface 
melting phenomena in atomic clusters should exhibit pro- 
nounced size effects, i. e. the melting behavior of two 
clusters, that differ in size only by one atom, can vary sig- 
nificantly. In particular, the (atomic) structure of closed 
(atom) shell clusters (magic clusters) is very sensitive to 
the addition of ad-atoms or vacancie a 16 ' 17 ' 18 ' 19 ' 20 . 

Our paper offers a systematic study of selected Al clus- 
ters in a range 12 < N < 562 near their melting transi- 
tion. Here, N is small enough, so that a simple extrapo- 
lation based on a the continuum theory is not applicable 
and new physics should emerge. Our most crucial ob- 
servation formulated in general terms: consider splitting 
the free energy of an TV-atom cluster into a bulk and a 
surface contribution 



F(T,N) = Nf B (T) + N 2 / 3 f s (T). 



(1) 
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Both terms, /b,Si depend on geometrical details of the 
cluster, i. e. we expect them to become strictly indepen- 
dent of N only in the limit A— >oo. Formally, /b.s are 
related to two reservoirs, called surface and bulk, with 
their own specific free energies. The reservoirs are cou- 
pled in the sense of the grand canonical ensemble, so 
they can exchange energy and particles. From this point 
of view, there is no reason why surface and bulk should 
have the same, or even a similar, melting temperature^ 1 - 
A reason why in metal clusters both temperatures tend 
to be strongly correlated with one another, nevertheless, 
is that the interatomic forces at the surface and inside 
the bulk are similar. 

Based on Monte-Carlo simulations employing semi- 
empirical Gupta potentials, we propose a general mecha- 
nism that can lead to a considerable splitting of the sur- 
face and bulk melting temperatures. Consider a closed 
shell cluster, e. g. AI13 or AI55, with icosahedral sym- 
metry. The outer shell of the II155 can accommodate an 
additional atom, an "impurity interstitial" , by replacing 
the fivefold ring structure surrounding an edge atom by 
a sixfold rosettelike ring. The formation of rosettelike 
structural excitations has been introduced already as a 
route to amorphisation of II155 systems^ 9 -, and we propose 
that is relevant for the binding of adatoms in II156 and as 
well. 

The impurity is mobile at the surface and its mo- 
tion strongly assists surface melting. This is, because 
the atoms inside the meandering deformation field are 
pushed away from their favorite high symmetry, low en- 
ergy site into a more shallow potential well at intermedi- 
ate position. Our explicit calculations strongly suggest, 
that by this mechanism the self diffusion of surface atoms 
can be dramatically enhanced indicating a significant re- 
duction in the activation barrier for diffusion and simi- 
larly also of the surface melting temperature. Since im- 
purities do not enter the next (second one counted from 
outside to inside) cluster shell, there is no correspond- 
ing reduction there, so that only two different melting 
transitions should be discriminated. 

This effect may have been seen in two recent exper- 
iments. Haberland et alM*Q have determined the la- 
tent heat and the melting entropy of sodium clusters, 
NaN, with A«50— 360. Their modeling of the data pro- 
vides an excellent phenomenological description assum- 
ing the premelting of the cluster surface for non-magical 
iV-values. The microscopic mechanism, that is responsi- 
ble for this lowering of the surfaces melting temperature, 
remained unspecified, however. In subsequent theoretical 
work the experimental melting temperatures have been 
reproduced quantitatively for a selected set of clusters^ 5 -. 
Moreover, the microscopies of premelting of sub-magic 
clusters has been already understood in terms of the dif- 
fusion of vacancie s 24 ' 26 . Still, the effect of ad-atoms has 
not been explicitly analyzed, and the relation to elasticity 
theory (strain) remained unexplored. 

Also, the specific heat of Aluminum cluster cations has 
been measured recently in the interval 49 < N < 63 by 



multicollision induced dissociation 27 . Interestingly, the 
specific heat data for a number of clusters shows signa- 
tures of multiple transitions, which have tentatively been 
interpreted as solid-liquid transitions at the surface that 
occur below the onset of melting^ 



II. METHOD 

To calculate the thermodynamics of an A-atom metal 
cluster, we employ a MC simulation in the canonical en- 
semble. Technical aspects of our procedure are described 
in detail in Ref. here, we focus on basic conceptual is- 
sues in order to provide the prerequisites necessary for a 
careful discussion of the numerical observables in Sees. IIIII 
and [TV] 

The potential energy of the metal cluster can efficiently 
be modeled by effective many-body potentials. We shall 
employ the Gupta potential (GP) 30 , which can be derived 
in the second moment approximation from a tight bind- 
ing mode l 31 ' 32 and which correctly describes the surface 
contraction observed in metals: 



V({r ij }) = J2 



N 

E 



Ae~ 




(2) 



Here, i and j are atom labels, Fy = ry/ro — 1, and 
Tij — \ri — Tj\ is the modulus of the distance between 
two atoms at positions and Tj. The parameters have 
been determined by fitting the experimental bulk lat- 
tice parameters and elastic moduli^ 3 , as A — 0.1221 eV, 
£ = 1.316 eV, p = 8.612, and q = 2.516 for Al. Dis- 
tances are measured in units of the bulk first neighbor 
distance ro = 2.864 A. A standard Metropolis algorithm 
is employe d 29 ' 34 ' 35 with boundary conditions imposed by 
a hard wall cube with linear dimension L: a shift of a sin- 
gle atom by a randomly chosen vector with a length taken 
from the interval [0, K,(T)ro] is offered with an associated 
change of the cluster energy AE. A temperature, T, is 
introduced via the probability, p, to accept such a step 
with p~ exp(- AE / k B T) . The parameter 0<k(T)<1 is 
chosen so that the acceptance rate is close to 50%; a typ- 
ical value at intermediate temperatures is k=0.25. The 
cluster is updated after each accepted move. Runs are 
performed with sampling rates of up to 8 x 10 steps per 
temperature and atom. 



Observables 

In order to characterize the thermodynamic state of 
the cluster we introduce the following observables: 

(i) The specific heat can be obtained from the ensemble 
averages of the potential energy V and its square, 



C_ 1 

~~ Afc 2 T 2 



((V 2 ) ~ (V) 2 ) + °- 



(3) 



3 



The kinetic contribution Cki n = 3/2fce per atom has been 
added here. (Note, that we treat the metal cluster as a 
gas of distinguishable particles rather than indistinguish- 
able ones. For the specific heat, a double counting prob- 
lem does not arise, since it is a second derivative of the 
free energy.) 

Invoking ergodicity, in practical calculations the en- 
semble average (. . .) is frequently combined with (or even 
replaced by) an average over the MC "time" , r, i.e. the 
total number of MC steps: 



r[K] 

464.2 928.4 



T[K] 
464.2 928.4 



r[K] 

464.2 928.4 
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(4) 



counting is started after equilibrating an initial configu- 
ration. The two body term O takes a value in the 
nth MC step. 

(ii) One may also introduce the rms pair index 



(5) 



to study the MC time evolution. The interest in this 
quantity stems from the following very general obser- 
vation: Consider a configuration space for an iV-body 
system which has the property, that several regions ex- 
ist where the free energy takes (local) minima. After a 
transient time interval ("warm up" or partial equilibra- 
tion period) the MC dynamics starts to explore one of 
these minima. There will be a typical time scale (corre- 
sponding to an activation energy) involved, after which 
the entire iV-body system migrates to a second, compet- 
ing minimum where the procedure repeats itself. Thus, 
the MC dynamics allows to study aspects of the energy 
landscape associated with the configurational space, like 
activation barriers. What has been described here for a 
general .ZV-body system remains equally valid for the pair 
of two particles, Eq. ([5]), embedded in an environment 
consisting of N— 2 other particles. 

In the limit r-^oo the evolution in MC "time" is er- 
godic, so the iV-atom cluster will explore all the phase 
space available. This implies, that at any temperature 
T > cluster atoms are deconfined: there is a finite time 
after which the i-atom has migrated from its initial posi- 
tion into any other given cluster site. Therefore, dij takes 
the same value, de for any given pair of atoms and is a 
unique function of the temperature and the particle num- 
ber: ds(T,N). Since all diffusion processes terminate at 
the cluster size independent of the diffusion constant (i. e. 
temperature) , g?b (T, N) has only a weak T-dependence 
about a mean value, that incorporates crude information 
about the overall cluster geometry (spherical vs quasi- 
oncdimensional), but nothing else. It is implied that the 
limit lie is usually not very sensitive to the melting tran- 
sition. 

For averages like the one defined in Eq. the ergod- 
icity theorem strictly holds only at r— >oo; at any finite r 
the value for dij(r) can vary between the different pairs 
of atoms If it so happens, that n different classes 
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FIG. 1: Characteristic temperatures for clusters of size N = 
12 (a)(d)(g), N = 13 (b)(e)(h), and N = 14 (c)(f)(i). Panels 
(a), (b), and (c) show the specific heat, panels (d), (e), and (f) 
the Berry parameter, while panels (g), (h), and (i) show the 
individual rms bond length fluctuations from Eq. {5]). Panels 
(a)-(f) show data for r=8 x 10 7 . In panels (g)-(i) dark symbols 
correspond to 8 x 10 7 MC steps while gray points are obtained 
with 4 x 10 7 MC steps per atom. 



of pairs exist, where each class just samples its own lo- 
cal minimum in phase space within r, then d%j (T) can 
develop n branches. 

Quite generally, in situations where different pairs sam- 
ple different sectors of phase space, the convergence with 
r can be increased by averaging over all the different 
pairs. In this spirit, we define one more average 



Sb(t) 



1 



N(N- 1) 



(6) 



which we will refer to as the "Berry parameter"—. It 
has a more rapid convergence behavior, 8b(t)— >ds, and 
therefore is easier to investigate in numerical simulations 
than the rms pair index. 



General properties of 8b and the rms pair index dij 

In order to illustrate the general properties of Berry 
parameter and pair index, we now consider as an exam- 
ple AI12, Ali 3 and Al i4 (Fig. [[]). We display the Berry 
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parameter Eq. ([6]) in Figs.QJd), (e), (f) for clusters with 
N = 12,13,14. 

At smallest temperatures only thermal vibrations of 
the atoms around a single site are observed within the 
MC window of time. Hence, the Berry parameter grows 
oc \/T reflecting the virial theorem applied to the har- 
monic oscillator. This asymptotic low temperature be- 
havior is clearly observed in the traces Fig. [TJ 

Also easily understood is the behavior at temperatures 
higher than Tg, at which the Berry parameter exhibits a 
very sharp jump. An estimate of 5b in this regime may 
be obtained, by taking the ground state geometry and 
calculating the average squared displacements by assum- 
ing ergodicity, i.e. that in the MC time evolution each 
occupation of allowed sites occurs with the same prob- 
ability. For example, in the case of AI13 one thus finds 
5b ~ 0.25 which agrees reasonably well with the data, 

Fig.ru 

The sharp increase of 5b at the intermediate temper- 
ature, Ts, signalizes the onset of cyclic, correlated ex- 
changes of (surface) atoms between their various posi- 
tions. (The weak irregularities, which are still visible, re- 
semble residual statistical noise.) In the spirit of our ear- 
lier discussion, we do not expect that T$ is independent 
of our observation time tM- In fact, the precise meaning 
of Tg is the following: at Ts the observation time r has 
been long enough, so that at T > Ts processes can be ob- 
served where atoms trade places with one another even 
though the probability p for this to happen may be ex- 
ponentially suppressed with a factor exp(— A/fceT). A 
denotes the corresponding activation energy which will 
in general exhibit a weak (i.e. non-singular) temperature 
dependency. 

Our argument shows that Ts itself can not immediately 
be identified with any intrinsic energy scale of the free 
cluster, like a surface melting temperature. The specific 
heat peaks only at a much higher temperature, Tc ^> Tg, 
which indicates the volume melting temperature of the 
cluster, see Figs. [TJa), (b), (c). 

It is possible to obtain an estimate of the activation en- 
ergy A from the way that Ts flows with the observation 
time t. Namely, one has A -1 ~ dTf 1 / dhifr) . Unfor- 
tunately, in order to obtain very accurate scaling with 
ln(r) the calculational effort goes well beyond what was 
achievable within this study. 

Further information is carried by the pair index dy, 
which is displayed in Figs.QJg), (h), (i) at two values of r. 
At N=13, a three branch structure is readily identified at 
T < Ts- The branches reflect the fact, that three kinds of 
atom pairs exist (center atom/shell atom, shell neighbors, 
shell next nearest neighbors), that have different distance 
fluctuations, see Fig. O Two branches are found with 
N=14 even at T > Ts- The upper one stems from the 
on shell pairings while the lower branch represents the 
mixed pairs, center atom/shell atom. Consistent with 
this picture, the lower branch contains 13 bonds, which 
are the 13 bonds between the center atom and the 13 
statistically equivalent surface atoms. This latter branch 



TABLE I: Nomenclature of relevant temperature scales of MC 
simulation with observation time r and the corresponding ac- 
tivation energies. 



T s 


onset of intralayer diffusion observed within r 


A 


To 


onset of interlayer diffusion observed within r 




Tc 


maximum of specific heat 





exists only in the intermediate temperature interval Ts < 
T < Tb>. The role of the temperature Td is similar to Ts, 
except that the associated activation energy, A; c , now 
is related to an exchange of inner atoms with the outer 
shell. At Td the Berry parameter, <5b, exhibits a second 
sharp increase. The discussion of the pair index allows us 
to attribute this increase as being due to the center atom 
now being deconfined within the MC window of time. 



III. SELF DIFFUSION 

The influence of closed shells on the cohesive energies 
of metal cluster o 7 ' 16 i 17 ' 18 i 29 in the gas phase and their 
melting point o 35 ' 38 has been a focus of research for quite 
some time. In this section we use the Berry parameter 
to investigate the onset of self diffusion of atoms within 
such systems. 



Ah 



The low temperature jump in the Berry parameter 
for AI14 in Fig. [T] has been identified previously as in- 
duced by the temporary absorption of the ad-atom on 
the AI13 icosahedral core structure into the surface, see 
Fig. [U 29 Here we witness the effect, that has been de- 
scribed in general terms already in the introduction. The 
ad-atom destabilizes the high symmetry surface of the 
AI13 cluster and therefore the activation barrier for self 
diffusion on the surface is reduced, by a factor of roughly 
T,s(Ali4)/7$(Ali3)~5 according to our calculation. 



4j j 




FIG. 2: Ground state configurations of AJ12 (left), AJ13 (cen- 
ter) and AJ14 (right) as obtained with Gupta-potentials. The 
hole has little impact on the position of the remaining atoms. 
By contrast, the ad-atom creates a docking site with a higher 
coordination number, which is incompatible with the three- 
fold symmetries of the docking sites of the unperturbed AI13 ■ 
Hence, AJ14 experiences significant strain. 
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FIG. 3: Potential energy statistics for Alu at ksT = 0.02 eV 
in the surface molten phase, (a) Convergence of the Berry pa- 
rameter as a function of MC steps, (b ) Sample of the potential 
energies of accepted configurations during the same MC run 
as in (a), (c) Histogram of the potential energy distribution 
function w(E) (2.4 x 10 6 energy values, 480 bins). 



In order to obtain a more detailed understanding of 
the thermodynamics of the system, in Fig. [TJi) the tem- 
perature dependence of the rms pair index is shown. In 
particular, we observe that Tu(A\ 1 4)>Ts(A\i S ). This 
suggests, that the core atom of AI14 must overcome a 
slightly increased barrier (as compared to AI13) to enter 
the (strained) outer shell. 

For illustration, the potential energy statistics for AI14 
are shown in Fig. [3] at an intermediate temperature 
(fceT = 0.02 eV), where the pair index exhibits two 
branches (Fig. Panel (a) shows the convergence be- 
havior of the Berry parameter as a function of r. Panel 

(b) displays potential energies of accepted configurations 
during the sample MC run. Panel (c) shows the potential 
energy distribution function, w(E). Consistent with our 
interpretation of the branching behavior and with earlier 
results for N114 the distribution shows no sign of phase 
coexistence 1 21 i 39 i 40 

The results presented here for AI14 are analogous 
to those for Pbi 4 with very similar characteristic 
temperatures^. They differ from Nii4^ and CU14— in 
so far, as for AI14 no ad-atom hopping is observed. 

The clusters with 15 < N < 18 show a behavior similar 
to AI14. Apparently, the presence of several additional 
ad-atoms has qualitatively a similar (destabilizing) effect 
on the surface as a single ad-atom. Panels (a), (b), and 

(c) of Fig.|4]show the individual bond length fluctuations 
Eq. ([5|) for Ali 5 _i7, respectively. For Ali 8 (not shown 
here, c.f. Figs. [7]) similar results are obtained. 
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FIG. 4: Temperature dependence of the individual bond fluc- 
tuations Eq. {5]) for clusters of size N — 15 (a), N = 16 (b), 
and N — 17 (c). The "solid" bonds with ckj < 0.1 correspond 
to those between the central atom and the surface atoms. 
In the surface molten phase all surface atoms are equivalent. 
Dark symbols correspond to 8 x 10 7 MC steps while gray open 
symbols are obtained with 2 x 10 7 MC steps per atom. 



atom, we now investigate the case Al 56 . 

Figure [H] shows the specific heat, the Berry parameter 
and the pair index for N = 55, 56, and 57, respectively. 
The dramatic suppression of T$ in Al 56 seen in <5b as com- 
pared to the closed shell case AI55 was already observed 
before^. 

To obtain more detailed information, we consider the 
pair index Fig. gji). The lower branch at T s < T < T D 
represents 331 pairs. They correspond to all pairings 
within the set of atoms that consists of the AI13 core to- 
gether with the 12 outer corner atoms of the (distorted) 
icosahedron - subset of confined atoms. The remaining 
31 surface atoms are deconfined. The upper branch rep- 
resents the cross-pairings between the two sets of atoms 
as well as the pairings between deconfined atoms. The 
data suggests, that edge atoms of surface facets are more 
mobile than the corner atoms. 

In order to understand the reason for this enhanced 
mobility, we display in Figure [5] the (Gupta-potential 
based) ground state structures of AI56 and AI55. The 
large full circles indicate the 12 corner atoms, which move 
very little at Tg < T < To. It is seen, that the ad-atom 




(b) 




B. Al 56 and AI57 

In order to see whether the lowering of activation bar- 
riers for surface diffusion may indeed be a typical phe- 
nomenon for closed shell configuration with one excess 



FIG. 5: Ground state configurations of AI55 (a) and Aha (b)as 
obtained with Gupta-potentials. Large full circles indicate the 
corner atoms, which remain "solid" in the partially surface 
molten state of AJ56. Small open circles (rosette structure^) 
show the broken fivefold symmetry caused be the 56 th atom 
absorbed into the surface of the icosahedron. 
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FIG. 6: Temperature dependence of the melting transition for 
clusters of size TV = 55 (a)(d)(g), N = 56 (b)(e)(h), and N = 
57 (c)(f)(i). Panels (a), (b), and (c) show the specific heat, 
panels (d), (e), and (f) the Berry parameter, while panels (g), 
(h), and (i) show the individual rms bond length fluctuations 
from Eq. {5]). Panels (a)-(d) show data for 8 x 10 7 MC steps 
per atom; in panels (g)-(i) dark symbols correspond to 8 x 10 7 
MC steps while gray open symbols are obtained with 4 x 10' 
MC steps per atom. 

integrates by promoting the environment of one of the 
corner atoms from its original 5 fold symmetry to an (ap- 
proximate) local 6 fold symmetry. This implies a sizable 
bond length mismatch of about l/5~20%, that creates 
a deformation field. Now, if an atom out of the six ring 
pushes away another atom to enter a new neighboring 
facet, then the entire deformation field will follow. Actu- 
ally, what really diffuses over the surface is the six-fold 
rosette structure (Fig. [5]). The energy barrier to be over- 
come in this self diffusion process is relatively low. This 
is because it is not necessary to first create a hole in the 
crystal lattice of the target facet for the rosette to move 
there. 

According to the scenario developed for AI56, two 
rosettes should decorate the surface of AI57. A naive 
expectation is, that these rosettes repel each another, 
because it is more difficult for a second rosette to enter 
an area that is under strain already from the presence of 
a first one. Since the two rosettes cover a large surface 
fraction of AI57, diffusion barriers should be increased 
again as compared to the "free" case AI56. Indeed, as 



can be seen in Fig. [6l the sharp increase observed for 
AI56 at very low temperatures in the Berry parameter 
almost disappears for AI57. It gives way to a very broad 
shoulder which is indicative of a large number of energy 
scales that is associated with a strongly disturbed ("dis- 
ordered") outer cluster shell. 

Experiments: In a recent experiment, a broad peak in 
the specific heat, C(T), of Al^" 6 and Alg" 7 has been mea- 
sured and analyzed i 27 ' 28 The authors were concluding, 
that Al^" 6 undergoes a separate transition even before the 
entire cluster starts to melt. The physical nature of the 
first transition could not be clearly resolved. It would 
be tempting to propose that our research supports the 
supposition of the authors that prcmelting of the sur- 
face is a possible candidate. Namely, since low diffusion 
barriers are usually also indicative of low melting tem- 
peratures, our scenario would suggest that premelting of 
the surface may occur well before melting of the bulk in 
Alg" 6 . However, such a direct application of our ideas to 
experiment a 27 ' 28 is not without difficulty and probably 
not indicated: there is a trend in the theoretical data, 
Fig. [fJl that the latent heat (integral under the peak in 
the specific heat) decreases from N=55 to N—57. This is 
plausible, because it takes less energy per atom to melt 
a structure under strain. By contrast, the experimen- 
tal trend is reverse and the latent heat increases for the 
series N=bb, 56, 57. The likely reason for this discrep- 
ancy is that Gupta-potentials cannot give a sufficiently 
realistic description of the thermodynamics of aluminum 
clusters, Aljy, (N=55, 56, 57) near the melting transition. 
That indeed difficulties exist even with the uncharged 
species, Al/v, is signalized by explicit density functional 
theory calculations, which show that the T=0 ground 
state conformation of AI55 is not icosahedra l 17 ' 19 contra- 
dicting what is found with Gupta potentials. This casts a 
doubt on the applicability of Gupta-potentials to the spe- 
cific metal aluminum for simulations of low temperature 
behavior. However, we would like to emphasize that our 
general findings remain valid for other materials with an 
isocahedral ground state as well as for those situations, 
where an approximate icosahedral symmetry is restored 
at slightly higher temperatures. 



IV. CHARACTERISTIC TEMPERATURES FOR 
ALat-CLUSTERS 

The delicate interplay between the special geometry of 
AI55 and the low activation barrier for ad-atom diffusion 
in AI56 and AI57 does not prompt the expectation, that 
this specific mechanism is ubiquitous in all its details. 
Nevertheless, there is a lesson to be learned about the 
more general case. One expects a lowering of the ac- 
tivation energy for surface diffusion of surface atoms of 
clusters with structures that derive from high symmetry 
parent states either by (i) punching in vacancies or (ii) 
by inserting ad-atoms into its outer shell. In such sys- 
tems, there is a possibility for atomic motions, which are 
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FIG. 7: Characteristic temperatures of AIn associated with self diffusion over the cluster size N; for the definitions ofT$,Tr>,Tc 
see Tab. U "a" labels a "confined" central atom, "d" a "confined" central dimer, and "t" a "confined" central trimer. (The 
onset of the evaporation transition in the investigated systems is found at larger temperatures r cvap > 0.12 eV^.j 



not just ring exchanges and which also do not require to 
break chemical bonds to first create a vacancy. Since low 
diffusion barriers are also indicative of low melting tem- 
peratures, here is a mechanism by which surface melting 
may become a process that should be distinguished, in 
principle, from the melting of the bulk. 

This analysis suggests, that the ratio Ts/Td tends to 
be large for structures with closed shells or subshells and 
much smaller otherwise. We have tested this idea by cal- 
culating Tq, the temperature at which the specific heat 
C(T) takes its maximum, T$ and Td (at fixed r) for a va- 
riety of different cluster sizes. Figure [7] comprises our re- 
sults which, we believe, support our general picture: the 
ratio oiTs/Tj) takes peak values at closed shell structures 
and much lower ones almost everywhere else (except for 
the smallest cluster sizes, where our previous analysis 
does not apply). 

We mention that the Al^v clusters with 14 < N < 18 
and N — 24 have maxima of the specific heat at tem- 
peratures larger than the Al bulk melting temperature 
of T bulk = 933 K = 0.0804 eV/k B . This observation is 
in line with the empirical investigation of the melting 
of small Sn 4 ^ and Ga^i clusters, which have revealed a 
possible stability of the solid phase of the particles be- 
yond the melting temperature of the bulk material. In 
these cases the high cluster melting temperature was in- 
terpreted as a consequence of the rigidity of the specific 
ground state structures of the clusters. This interpre- 
tation found support from microcanonical molecular dy- 
namics (MD) calculations for C, Si, Ge, and Sn clusters 45 
as well as for isokinetic MD investigations of Sn m 27 i 46 
and Gai3, Gai^ particles. 

We briefly touch upon the limit of large clusters, 
N > 100. There, we observe that T D —>Tc, while T s 



does not appear to follow this trend. To understand this 
behavior, recall that Td is the temperature at which our 
MC time has become long enough, so we can observe 
an exchange of particles between the outer shell and its 
inner neighbor. Then, Tb«Tc implies that intershell ex- 
change cannot be observed - even with our very long 
observation times - unless we actually heat up the entire 
cluster to melt. This behavior is consistent with what 
one would expect for a macroscopic single crystal grain: 
as long as the crystalline structure of the surface is intact 
(T <C T me it) one has A <C Ai G and interlayer-diffusion is 
strongly suppressed. However, after melting the surface 
layer no longer forces the atoms of its neighboring inner 
shell into a crystalline structure. Hence, intra self diffu- 
sion within the second layer becomes (almost) as cheap 
energetically as was diffusion in the first layer, before. 
Therefore, the second layer melts immediately after the 
first one, so at melting one has for the "effective" acti- 
vation energies: A « Ai Q ~ T mo i t . In other words: we 
recover the standard "continuous melting" scenario. 

V. CONCLUSION 

In this paper we have argued, that certain metal clus- 
ters with an intermediate size may exhibit a property 
that cannot be found in bulk materials: the activation 
energies for surface diffusion, A, and interlayer diffu- 
sion, Aj , are substantially different from one another: 
Ai ^> A. Since these activation energies are closely tied 
to melting temperatures, one expects that the outermost 
surface layer can exhibit its own melting transition, which 
is well separated from the bulk. The continuous melting 
of the cluster core should start only at much higher tern- 
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peratures. 

The most dramatic decrease of A has been found with 
AI56. In this case, it is the "frustration" of atomic bonds 
that originates from implanting an ad-atom into a closed 
shell system, which produces the effect. A related mech- 
anism leads to a decrease in the surface melting tem- 
perature also for AI14 and is expected to be active in 
clusters, where N is slightly above some magical (closed 
atom shell) value TV*. We mention, that a distortion of 
the outer cluster shell is also present in clusters with N 
slightly below N* . However, a hole is usually accom- 
modated more easily than an ad-atom and therefore the 
decrease tends to be asymmetric: it is typically stronger 
for N > N* as compared to N < N*. 

Monatomic, macroscopic and planar metal surfaces, 
that face the vacuum, do not easily allow for the frus- 
tration of surface bonds that we have observed with the 
metal clusters. An attempt to locally implant concen- 
tration of ad-atoms into macroscopic surfaces beyond a 
certain threshold would result in a metastable state that 
eventually would transform into another state without 
frustration, where the ad-atoms would have undergone 
island formation. 

On the other hand, our research suggests the design 
of materials with a surface melting temperature that is 



strongly diminished and separates from the bulk melt- 
ing temperature by a controllable amount. The idea is 
to employ a monatomic core and a biatomic shell struc- 
ture. The purpose of pressing foreign atoms into the 
outermost layer of the host material's crystal is to create 
local strain fields. Since strain reduces the local melt- 
ing temperature, heating up such a system could create 
puddles of molten host material on top of the solid, bulk 
core. Clearly, the combination of host and implantation 
materials should satisfy at least two conditions: (i) the 
implantations should have a high solubility in the host 
material, but (ii) they should not easily diffuse away from 
the surface into the bulk of the crystal, either. Whether 
indeed a combination of suitable materials can be found, 
this we have to leave for future research. 
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